Dimerized phase and transitions in a spatially anisotropic square lattice 

ant iferromagnet 



Oleg A. Starykh^ and Leon Balents^ 
^Department of Physics and Astronomy, Hofstra University, Hempstead, NY 11549 
^Department of Physics, University of California, Santa Barbara, CA 93106-4030 



o 
o 

(N 
m 



We investigate the spatially anisotropic square lattice quantum antiferromagnet. The model 
describes isotropic spin- 1/2 Heisenberg chains (exchange constant J) coupled antiferromagnetically 
in the transverse (Jx) and diagonal (Jx), with respect to the chain, directions. Classically, the 
model admits two ordered ground states - with antiferromagnetic and ferromagnetic inter-chain 
spin correlations - separated by a first order phase transition at J± = 2 Jx . We show that in 
the quantum model this transition splits into two, revealing an intermediate quantum-disordered 
columnar dimer phase, both in two dimensions and in a simpler two-leg ladder version. We describe 
quantum-critical points separating this spontaneously dimerized phase from classical ones. 
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An interplay between geometric frustration and quan- 
tum fluctuations is at the heart of intensive current in- 
vestigations into the nature of possible S'[/(2)-invariant 
Mott insulators and quantum phase transitions. The ab- 
sence of any "natural" small parameter, however, makes 
a traditional perturbative approach difficult. Aside from 
direct numerical attacks involving exact diagonalization 
and quantum Monte-Carlo techniques, this feature forces 
one to extend the parameter space of frustrated mag- 
netic models and explore various "corners" of the result- 
ing phase diagram in the hope of gaining new insights into 
the physically relevant region of system's parameters. In 
this paper, we explore a simple spatially anisotropic frus- 
trated spin model which allows the application of power- 
ful analytical methods borrowed from one dimension. 

The model we consider - the spatially anisotropic 
square lattice antiferromagnet - is a quasi-one- 
dimensional generalization of the well-studied square lat- 
tice antiferromagnet with frustrating antiferromagnetic 
exchange along diagonals, also known as the Ji — J2 
model. It describes a collection of antiferromagnetic 
spin— chains with exchange constant J running along 
horizontal (chain) direction. The chains are interacting 
with their nearest neighbors via weak antiferromagnetic 
spin exchange in the transverse {J±) and diagonal (Jx), 
with respect to the chain, directions. Recently this model 
was investigated by Nersesyan and Tsvelik (NT) (ll| who 
predicted a novel RVB-like phase with deconfined mas- 
sive spinous at the special ratio of microscopic exchange 
constants: J± = 2Jx <C J. We shall comment on this 
work and subsequent numerical investigations 0, 
the end of our paper. Let integer index n numerate sites 
along the chain while index m numerates chains. Then 
the Hamiltonian is H = J2m + where Hm^ is the 
standard Heisenberg Hamiltonian of the m-th chain and 
the interchain interaction reads 
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summation is implied. The classical {S = 00) phase dia- 
gram is simple: for J± > 2 Jx ordering of the spins in the 
transverse to chain direction is antiferromagnetic (AFM 
phase), whereas for J± < 2Jx it is ferromagnetic (FM 
phase); spins order antiferromagnetically along chains. 
In the isotropic Ji — J2 model these orderings correspond 
to the Neel and four-sublattice phases, respectively. The 
line J_L = 2 Jx describes a first order transition between 
these two phases. Along this line the interchain interac- 
tion can be written in a more suggestive form 
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Our goal is to understand the phase diagram of the quan- 
tum = 1/2 model in the limit of weak interchain cou- 
plings J±, Jx J- Taking the continuum limit along 
the chain direction (x — na, a is the lattice spacing), the 
spin operator is decomposed as 



+ ^™,L(^) + (-ir^™(^) (3) 



in terms of chiral components J^ of the uniform spin 
magnetization of the m-th chain, and the staggered spin 
magnetization iV°j (scaling dimension 1/2). The uniform 
spin magnetization Jm = ^m.L the conserved 

spin current (scaling dimension 1) and a generator of 
0(3) spin rotations. The two fields in ^ are connected 
via the following operator product expansion (OPE) 
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Here a = x^y,z is the vector index over which implicit 
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Here the upper/lower signs on the right hand side ap- 
ply for the right/left (R/L) moving currents, respectively, 
and V = TT Ja/2 is the spin velocity. Of major importance 
is the appearance of the staggered dimerization operator 
e (scaling dimension 1/2). It is the continumn limit of 
the lattice operator involving the scalar product of two 
neighboring spins, e™(x) ~ (-l)"5'^,™S','J+i,™- 

The spin algebra encoded in Q and well-known OPE 
of the chiral components of the spin current 2] follow 
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from the fact that the Hamiltonian density of the isolated 
Heisenberg chain is described by the SUi{2) WZW model 
perturbed by a marginally irrelevant backscattering term 
(.9bs = 0(J)) 
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The c-functions appearing in the OPE above are just 
the Green's functions of right- and left-moving fermions 
which are governed by the Hamiltonian ^ with gbs = 0. 

Having exposed the structure of the unperturbed the- 
ory, we now turn to the inter-chain interaction JQ). Its 
continuum low-energy Hamiltonian density contains all 
terms allowed by the symmetries of the lattice model 
(these include reflection with respect to transverse direc- 
tion and translation by one lattice constant) 
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Here dx denotes derivative with respect to x and the bare 
coupling constants follow from substituting (p!]) in l|T|) 

91 = {J± - 2Jx)a,.g2 = (J± + 2Jx)a,g3 = -|-'54 = 0. 

(7) 

The scaling dimensions of the first, second, and third in- 
teraction terms in © are 1,2, and 3, respectively. The 
last term (scaling dimension 1) is included here because 
it respects the symmetries of the lattice model . Al- 
though its bare coupling constant is zero, quantum fluc- 
tuations can (and will) generate some finite 34. 

At this level the phase diagrams of the quantum model, 
© with 34 = 0, and classical one, ^ with S = 00, are 
identical. As long as J± ^ 2 Jx , © is dominated by the 
strongly relevant interaction of staggered magnetizations 
which drives the system into one of the classically ordered 
phases: AFM for gi > and FM for gi < 0. 

However, the transitional region J± « 2Jx requires 
a closer look because there the amplitudes of the both 
relevant terms are zero. To begin, let us fine-tune to the 
point gi = 0. On the lattice this corresponds to ||2J), 
whose continuum limit is given by (|SJ) with gi = 94 = 0, 

Vx = ^ 92J^J:n+i + g:ia^dxN^dxN'^^^. (8) 
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The situation is now apparently controlled by the 
marginally relevant (52 > 0) current-current interaction 
which seems to suggest that the strongly irrelevant g^- 
term does not play any role. Recall that the (Kac- 
Moody) algebra of spin currents ^^/^ is closed in the 
sense that OPE of two like currents produces another 
current of the same chirality . Thus Vx with gy, ~ 
will "reproduce" itself in every order of the expansion in 
powers of the coupling constant 52- This peculiar behav- 
ior is destroyed by 33 7^ (or any other equally or even 
more irrelevant term). To see so, we expand the partition 



function to second order in Vx . As follows from the dis- 
cussion above, one should concentrate on the cross-term 
(« .9253) 
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Applying the OPE I^J to the same-chain operators at 
nearby points, differentiating the c-function pre-factors 
with respect to x and performing elementary integration 
over the relative coordinates, one obtains a relevant cor- 
rection to the Hamiltonian of the same form as The 
coefficients in ((TJ are thereby replaced by 
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51 = [Ji- - 2Jx + -^)a, 
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while two other couplings (72,3 remain unchanged. Ex- 
actly the same result (O can be obtained by perform- 
ing the standard renormalization group (RG) analysis of 
(O and ®. There one obtains five coupled non-linear 
differential equations for 51-4 and ghs- In the transi- 
tional region, corresponding to 31,4 = 0{J^/J) <^ J^, 
RG equations decouple and admit a simple analytical so- 
lution, which reproduces ©, for £ of the order 1 (here i 
parametrizes the short-distance cutoff: = aoe^). Our 
derivation of the renormalized coupling constants is very 
similar in spirit to the recent calculation of second-order 
corrections to coupling constants of two relevant and 
competing interactions in the extended Hubbard model 

i- 

What are the phases of V with relevant couplings given 
by lO? We begin our analysis by studying two cou- 
pled chains, that is, frustrated spin ladder (m — 1,2 
in Q). Hamiltonian density of the ladder Tiiaddcr = 
Ti-f^ + h'^^ + V can be recast as a theory of four massive 
real (Majorana) fermions, see Ch. 21 of 0, 
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where JTimarg, originating from gbs m (|SJl and 52 in 0, 
describes residual marginal interactions between Majo- 
rana fermions 
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with g± — g2 ± ghs- Ignoring SHmmg, the first three 
Majorana fermions form a triplet with the mass 



whereas the fourth fcrmion ■C^/^:, is a singlet with 



(12) 



(3.91 +94) 3A2 _^Jl. .... 

TOs = = —{JL-2Jy,+^T-). (13) 



27ra 



27r 



3 



Both masses receive finite logarithmic corrections from 



which we omit in the following. Here A ^ 1 is 



the expectation value of the charge operator which mul- 
tiplies abelian bosonization expressions for the staggered 
magnetization and dimcrization operators Q 

N„i = (cos \/2TT9m , sin \/2TT9m , — sin \/2TT(p,n ) , 

TTfl 

Cm = —C0sV2TTipm- (14) 

na 

The Hamiltonian (|10|) predicts two transitions fs] which 
happen when one of the renormalized masses becomes 
zero. 

(?) mt = 0: at this point triplet excitations become gap- 
less. Integrating out massive fermion we find that 
^ladder is the Hamiltouian of the SU2{2) WZW model 
with central charge c = 3/2. The nature of the transition 
is determined by the sign of remaining marginal coupling 
92 — 9hs- Its positive initial value implies marginally rel- 
evant flow to the strong coupling, and hence, first order 
phase transition. For the negative initial value one ob- 
tains marginally irrelevant flow of the coupling constant 
to zero, and hence, second order transition. 
(a) nis — 0: singlet excitations become massless. This is 
a continuous Z2 (Ising) transition. 

Properties of the various phases are most conveniently 
understood from the bosonized form of the relevant part 
of interchain interaction lO, Vrci, which corresponds to 
mass terms in ()10|l . In terms of symmetric and anti- 
symmetric combinations of bosonic fields on two chains, 
f± = {'Pi ^ ^2)1^/^ (and similarly for B±) it reads 



Vrcl — 



A2 



2gi cos V47r0_ + (gi + 174) cos 

(15) 



2(7ra)2 

-(31 -54)cos\/47r(p 



Depending on the ratio gi/gi four phases are possible: 

(a) gi — > —00 is the Haldane (effective spin-1 chain) 
phase; mt < 0, > and = \/7r/2, ^- = 0. 
Inter-chain spin correlations are ferromagnetic. This is 
the ladder analog of the classical FM phase. 

(b) gi +00 is the rung-singlet phase; mt > 0, m^ < 
and Lp+ — 0, 9- = V^/"^- Here inter-chain spin corre- 
lations are antiferromagnetic, which corresponds to the 
classical AFM phase. 

(c) gi —I- —00 is the columnar dimcr phase (DC); m^^t > 
and (^-|_ — 0, = 0. 

(d) gi +00 is the staggered dimer phase (DS); ms^t < 
and if^ = if- = 

The last two phases are spontaneously dimerized and 
do not have classical prototypes. Using (|14|l we find 
that both dimerized phases break translational symme- 
try (along the chain direction) but their ordering pat- 
terns are different. In the DC phase (ei) = (£2), which 
describes columnar dimer long-range order. In the DS 
phase (ei) — —(£2), which makes it the staggered dimer 
state. Observe that microscopic couplings 10 choose 
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FIG. 1: Phase diagram of the frustrated ladder. Line AB 
marks the parameters corresponding to the lattice model. 



columnar dimer as the intermediate phase between the 
Haldane and rung-singlet phases of the frustrated ladder. 

This feature was missed by (fil, 1^ who described J± = 
2Jx point of the lattice ladder by © with g^ = 0. Our 
analysis shows that multicritical point where chains are 
coupled only by the marginal current-current interaction 
requires additional fine-tuning of both gi , 514 (equally, of 
ms,mt) to zero. A more detailed understanding can 
be gained by considering the model ^2 > and expo- 
nentially small relevant couplings gi^ (more specifically 
1/1 In gi,4| ^ 52) - though this is not realized by the lat- 
tice model. In this case the NT analysis 1] can be taken 
over, treating gi^^ as small perturbations. A scmiclassi- 
cal analysis for gi^ = shows that the NT model has 
a manifold of ground states in which minima of type (a) 
and (b) above - i.e. Haldane and rung-singlet states - 
are degenerate. Perturbing with H15|) above splits the en- 
ergy densities of these two degenerate states, so that the 
gi — 3gi becomes a first order transition line (dashed line 
through the origin in Fig. Furthermore, the spinons 
of the NT phase correspond to domain walls between 
these two coexisting states. Remembering that the finite 
length spin chains in the Haldane phase are characterized 
by spin-1/2 end states, the spinons become physically 
transparent as mobile versions of these. Away from the 
first order line, such spinons are therefore conGncd (i.e. 
cost an energy linear in system size), due to the energy 
cost to create a domain of the disfavored phase. Thus 
although it does not occur in the lattice model studied 
here, the NT state and its spinons can exist as a first 
order line in a generalized ladder model. 

Having understood the ladder, we tackle the full two- 
dimensional problem. We treat the relevant inter-chain 
interactions giA (taking 172.3 = 0) using a self-consistent 
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mean- field approximation!^ |^. Assuming symmetry 
breaking along z-dircction, 

V = -z(|<?i|(iV,l)7V^, + |<?4|(e™)e™) (16) 

where we used {N^^j^^) — ±{N^) for gi negative (pos- 
itive), and z = 2 is the number of neighboring chains. 

Using (fnji and bosonized form of Tim^ (with ^bs — 0), 
rescaling euclidian time t — y/v and introducing dimen- 
sionless euclidian distance r = (x,y)/a, we arrive at the 
following single-chain sine-Gordon action 




where dimensionless parameters 

p^f^(sinV2^^), g^f^(cosV2^^) (18) 

are effective staggered magnetization and dimerization 
fields, respectively. These averages are found by differen- 
tiating the free energy density F{p, q) = — \nZ/V, Z — 
Tr(e^^'°') with respect to p and q, respectively. Clearly 
from IjlTI) , the mean-field free energy can depend on p, q 
only through k = \J p^ + q^ . One can therefore take 
P{V-, q) — f^) and take advantage of exact results 
for the standard sine-Gordon action[lOj. In particular 
dF/dn = — ciK^/^, where the numerical constant ci fol- 
lows from equations (10-14) of ref.p!o| 

_ tan(7r/6) / 2r(l/6) / 7rr(3/4) x4/3 
' 3 Vv^r(2/3)^ V2r(l/4)^ • ^ ' 

Thus our self-consistent equations 

(sin \/27r(p) = , {cosV2TTip) = -^^^ (20) 

become simple algebraic ones. They predict two tran- 
sitions, at gi = ±|(74|, where both order parameters 
(N^,em PtQ) are non-zero. This makes them first- 
order transitions - the first order nature being likely an 
artifact of the mean-field approximation. These two tran- 
sitions separate three phases, indicated in Fig. |21 For 



\9i\ > I54I: P 7^ 0, g = 0, so the system is AFM or FM, 
depending on the sign of 171. For \gi\ < \g4\, one finds 
the columnar dimer phase (rather than staggered since 
g4 < 0). Note that obtained phase diagram (Fig. O is 
in qualitative agreement with that of the ladder, which 
provides a posteriori support of our mean-field treatment. 
Moreover, it agrees very well with the phase diagram of 
the lattice model obtained in the exact diagonaliza- 
tion study 0. We should add that DC phase was not 
observed in DMRG study ^ , presumably due to strong 
finite-size effects at small J±/J. Now, since exactly this 
type of dimer ordering is known to take place in the spa- 
tially isotropic Ji — J2 model 0, it is natural to con- 
clude that the DC phase found here extends all way up 
to J± = J. 

Jx/J 




FIG. 2: Phase diagram of the two-dimensional model. Solid 
(dashed) hnes denote phase boundaries of the quantum (clas- 
sical) model. 
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